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■ Abstract 

\ It is shown that the existing data on two-body B decays, some of them only 

' upper Umits, are precise enough to perform an isospin analysis to extract the phase 

. shifts due to final state interaction. Unlike charm decays, no significant final state 
interaction is observed in decays B Dtt, Dp, and Z)*7r supporting the factoriza- 

Q^' tion hypothesis in these decays. From the isospin amplitudes obtained, we extract 

■ the ratio a^jax, where a\ and 02 are the coefficients in the factorized effective Hamil- 
D \ tonian. 

> 

^ ■ The idea of factorization has been used for evaluating non-leptonic weak decays ever 

■ since Schwinger employed it to show that the A/ = 3/2 kaon decay rate is consistent with 
the corresponding semileptonic decay rate JI], 0. Even though the method was originally 
thought to be useful only for order of magnitude estimations, it has been extensively 
applied to heavy hadron decays 0, |, |^, ^ with results varying from mixed to reasonable. 
Since factorization is thus far virtually the only way to quantitatively calculate exclusive 
non-leptonic rates, it is important to test the validity of the hypothesis whenever possible; 
in fact, it has been a subject of a number of studies 0. 

We take the decay B D^ti~ as an example which can occur by the 4-fermion 
operator p 

a^{du),,{chY + a2{cu)^{dhY (1) 

where the notation is ijlq')^ = ^7^(1 — 75)?' and ai^2 are real coefficients which are related 
to the Wilson coefficients Ci and C2 by 

ai = Ci + eC2 a2 = C2 + eCi. (2) 
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The parameter ^ is sometimes called the color suppression factor and naively expected to 
be 1/3. At the mass scale of 6-quark, the leading- logarithm approximation gives 

Ci{mb) = 1.11 Ciinib) = -0.26. (3) 

The basic idea of factorization is that the pion is generated from vacuum by the 
current operator {du)^^ and the transition B ^ D is caused by the current operator (cb)^, 
and that they occur independently. In terms of matrix element, it amounts to the fact 
that it can be written in a factorized form: 

Amp(B' D^n-) = ^V:,V^ba,{n-\(dum{D^\{-cb),\B'). (4) 

Such assumption of factorization has been shown to be correct to the zeroth order in 
1/A^ expansion 0. Also, it has been argued intuitively that factorization should hold for 
energetic tow-body decays p, |12| based on two observations: 1) when the ud pair escapes 
the color field around the b quark, it is highly energetic (~ 2.5 GeV) and, since it has to 
eventually form a pion, the pair is collinear and close together with the total color being 
zero. Thus, we expect that the pair will escape the color field without much interaction 
('color transparency'). 2) By the time the pion is formed, it will be well outside the 
color field (again due to the high energy of the pion); thus, there will be little final state 
interaction (FSI) between the D meson and the pion. The same arguments had been 
used for p and ijj productions in hard scatterings [|13], and recently put forward further 



by Dugan and Grinstein in the framework of heavy quark effective theory [|14|. 

Since factorization assumes independence of the pion formation and the B to D tran- 
sition, FSI between them is antithesis to factorization. Often, however, factorization and 
FSI are combined to be compared with data assuming that the factorization calculation 
correctly estimates the amplitude 'just before' FSI takes place P, |l^. This procedure has 
worked reasonably well for charm decays, and we will effectively employ it later when we 
extract 02/01. It is worth keeping in mind, however, that such treatment is not justified 
in the 1/A^ expansion ^ and also it is not well defined exactly where factorization ends 
and FSI begins. 

In the isospin analysis of the charm decays D Ktc it is found that the FSI 
phase shift between isospin 1/2 amplitude and isospin 3/2 amplitude is large (~ 77°). 
Furthermore, after the effect of FSI is removed as described above, the coefficients oi and 
02 are found to be nearly identical to the values of Wilson coefficients Ci and C2 evaluated 
at the charm mass scale, corresponding to ^ ~ 0. This has led to the rule of so-called 
'discarding 1/N terms' H and prompted further theoretical studies based on QCD sum 



rules [|rT|. A recent analysis [10|, however, indicates that the situation is quite different 



for B decays giving the value of ^ around 1/2 to 1/3. In this study, we will see that the 
existing data leads us to conclude that the FSI phase shifts are small for the B decays 
B — s> Dtt,Dp, and D*tt, and that the value of ^ remains to be l/2~l/3 even after the 
effect of FSI is taken out. 
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The Hamiltonian responsible for B Dtt decays has isospin |1, — 1), and this leads 
to the following isospin relations: 



where A+- = Amp{B^ D+n-), = Amp(B^ ^ D'^n^), A^' = Amp{B- D 
and A3/2, v4i/2 are the isospin 3/2 and 1/2 amplitudes, respectively. There are three 
unknown parameters: IA3/2I, |^i/2|, and 5 = aig{A3/2/Ai/2). Since there are three mea- 
surements of decay rates (namely, |y4+_p, |v4ooP, and IAq-P), one can solve for the three 
unknowns: 

lAsl = 



TT 



I An 



l^il = \/|A+-P + |AooP-^|^ (6) 

. |A+_p + |Aoop-|Ao_|V3 

cos = I = 

^||Ao-|2(|A+_|2+|Aoo|2-|Ao-|V3) 

Note also that the isospin relations (||) can be expressed as a single triangle relation: 

A+- + v^A°° = (7) 

The same relations hold for the decays B D*tt and B Dp. For the decay 
B D*p the same relations hold separately for each helicity amplitude. Since there is no 
a priori reason to believe that the polarization is the same for D*'^p~' , D*^p^, and D*^p'^, 
and since there is not enough data to separate the helicity amplitudes, we will not include 
D*p mode in this analysis. 

Amplitudes calculated by factorization naturally satisfy the triangle isospin relation 

(^. This can be seen from the expression (^) and corresponding factorized forms for A^ 

and Aqq: 

^ ^K;:iK(,«2{C°l(a<)''|0)(7r»|(3(,)„|B°) 

V2 



°- - a,{n-\(udr\0){D'\(cb),\B-) 



+a2{D"\icur\0){n-\idb)^\B-) 



and noting that (from isospin symmetry) 

{D^cbUB-] 
V2{n'\(db),\B') = {n-\(db)^\B-). 



{D'm,\B-) = {D^\(cb),\B') 
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In fact, the isospin amplitudes are explicitly given by 

Al = ^lC,V;t^(Afi + A/2) 



or 



where 



^3/2 1 + M2/M1 

2-M2/M1 ^ ^ 

Ml = ai(7r-|(rfu)^|0)(D+|(c6)^|#) (10) 
M2 = a2(/^°|(cM),JO)(7r-|(56)^|fi-) 

As before, the same relations (|8|-|T0D hold for decays D*7r, Dp, and each helicity state of 
D*p. 

If factorization is assumed, the three amplitudes ,^005 and Aq_ are relatively 

real, and thus the triangle reduces to a hne. Therefore, if the measured decay rates 
are exactly as expected from factorization as prescribed above, then the isospin analysis 
is guaranteed to give 5 = 0. Actual measurements, however, are always associated with 
errors, and the main point of this article is in showing that meaningful isospin analyses 
can be performed even though only upper limits are available for some of the decay modes. 

The isospin amphtudes cannot be uniquely given by factorization since it depends 
on decay constants and form factors through (|To|) . If we use the model of Bauer, Stech 



and Wirbel g [T^] together with = 132 MeV and fo = 220 MeV, we obtain 



M2 

Ml 



1.23 02/01 (Dn) 

1.30 02/01 iD*TT) (11) 
O.6602/01 (Dp) 



where the 01,02 are the coefficients appearing in the effective Hamiltonian (|l|). This can 
be substituted in (^ to obtain the expected ratio of isospin amplitudes, or if the ratio 
^3/2/^1/2 is known, 02/01 can be extracted from 

M2 _ ^ A3/2/A1/2 - I/V2 
Ml A3/2/A1/2 + V2 ' 

The table 1. shows the current available measurements for the relevant decay modes 



| T0|| . We will take the statistical errors only, and assume that the life times of B and B 



are the same: 

r(S°) = r{B-) = l.lSps. (13) 

The upper limits are converted to a gaussian distribution centered at zero by setting the 
r.m.s. of the gaussian to (upper limit)/1.64. 
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mode (%) 


fi" mode (%) 


B- mode (%) 


D+Ti' 0.29 ±0.04 
±0.03 ±0.05 

D*+Ti- 0.26 ±0.03 
±0.03 ±0.01 

D+p- 0.81 ±0.11 
±0.12 ±0.13 


D^Ti^ < 0.035 
L)*°7r° < 0.072 
< 0.042 


D^TT- 0.55 ±0.04 
±0.03 ±0.02 

D*^Ti- 0.49 ±0.07 
±0.06 ±0.03 

D^p- 1.35 ±0.12 
±0.12 ±0.04 



Table 1: Branching ratios measured by CLEO. The first error is statistical, the second 
error is systematic, and the third error is due to uncertainties in D branching ratios. The 
upper hmits are 90% confidence levels. 



Table |^ shows the solution for IA3/2I, |^i/2| and cos 5 using the formulae (|^). For 
each mode, cos 5 is consistent with unity indicating that there is no phase shifts due 
to final state interaction. In this analytical method, however, the range of cos 5 is not 
constrained to within ±1. In order to take the constraint into account properly, we will 
use the maximum likelihood method. The likelihood function for |v43/2|, |^i/2|5 and cos 5 
is given by 

i=+-,oo,o- v27ro-i V ) 

where Fj and cxj are the measured decay rate and its error, F^ is the decay rate calculated 
from the amplitude given by according to the standard formula 

T.0_ V ,^i|2^ ^^^^ 



and the normalization factor is added to make the integral over the allowed region of 
1^3/2!, 1^1/2!, and cos (5 to be unity. Result of the fit is also shown in Table ||. Because 
of the constraint on cos 5, the errors are generally better than those of the analytical 
solutions. For each mode, the most likely value for cos 5 was unity. Figure |T] show the 1,2 
and 3 sigma contours for IA3/2I vs IA1/2I and IA3/2I vs cos 5. It is seen that the parameters 
are not strongly correlated. 

Also shown in Table || are the ratio of isospin amplitudes IA3/2M1/2I and 02/01 
extracted using equations (^T|-[T^). For each mode, the ratio 02/01 is positive which is a 
consequence of IA1/2I < v^|^3/2|- Averaging over the three modes, we obtain 

— = ±0.25 ±0.05 (16) 

Oi 

or from (||) the corresponding color suppression factor ^ is 

^ = 0.45 ±0.04 (17) 
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Figure 1: One, two and three sigma contours for the 2-dimentional plots of the hkehhood 
function: IA3/2I vs |^i/2| (a) and |>l3/2| vs cos 5 (b). 
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Dn 


D*7l 


Do 


analytical 


^3/2 


(10-5 

(10-5 


Vjfc! V J 


556 ± 021 


533 ± 038 


886 ± 040 


solution 


^/2| 


GeV) 


0.425 ± 0.095 


0.411 ± 0.127 


0.792 ± 0.133 




COS 6 




1.20 ±0.28 


1.19 ±0.61 


1.11 ±0.13 


maximum 


^3/2 


(10-5 
(10-5 


GeV) 


0.550 ±0.020 


0.527 ±0.037 


0.862 ±0.037 


likelihood 


^1/2! 


GeV) 


0.503 ±0.057 


0.494 ±0.065 


0.907 ±0.084 






COS 5* 




> 0.82 


> 0.57 


> 0.92 


^43/2Ml/2 




1.09 ±0.13 


1.07 ±0.16 


0.95 ±0.10 


02/01 






0.25 ±0.07 


0.22 ±0.08 


0.31 ±0.11 



* In all cases the most likely value for cos 6 is unity. The lower limits are at 90% confidence 
level. 

Table 2: Analytical solutions and results of the maximum likelihood fit for the isospin am- 
plitudes and their relative phase angle. Also given are the ratio of the isospin amplitudes 
and 02/01 derived therefrom (using the result of the maximum likelihood fit). 



which is consistent with the analysis of Ref |]T0| where the decay rates were fit to the 
model by Bauer, Stech and Wirbel without taking out the final state interaction. This, 
however, cannot be considered to be an independent confirmation of the positive value of 
02/01 since the two analyses are highly correlated. 

In summary, we have performed an isospin analysis on two-body B decays and found 
that the phase shifts by final state interaction are small in stark contrast to the case of 
charm decays. By fitting to the obtained isospin amplitudes, we have also seen that the 
effect of removing the final state interaction does not alter the observation that the ratio 
O2/01 is positive. 
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